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ABSTRACT 

We present observations of Arp 147, a galaxy system comprising a collisionally- 
created ring galaxy and an early-type galaxy, using the Oxford SWIFT integral field 
spectrograph (IFS) at the 200-inch Hale telescope. We derive spatially resolved kine- 
matics from the IFS data and use these to study the interaction between the two 
galaxies. We find the edge-to-edge expansion velocity of the ring is 225 ± 8kms _1 , 
implying an upper limit on the timescale for the collision of 50 Myrs. We also calcu- 
late that the angle of impact for the collision is between 33° — 54°, where 0° would 
imply a perpendicular collision. The ring galaxy is strongly star-forming with the star 
formation likely to have been triggered by the collision between the two galaxies. 

We also measure some key physical parameters in an integrated and spatially 
resolved manner for the ring galaxy. Using the observed B— I colours and the Ha 
equivalent widths, we conclude that two stellar components (a young and an old 
population) are required everywhere in the ring to simultaneously match both observed 
quantities. We are able to constrain the age range, light and mass fractions of the young 
star formation activity in the ring, finding a modest age range, a light fraction of less 
than a third, and a negligible (<1%) mass fraction. We postulate that the redder 
colours observed in the south-east corner of the ring galaxy could correspond to the 
nuclear bulge of the original disk galaxy from which the ring was created, consistent 
with the stellar mass in the south-east quadrant being 30-50% of the total. The ring 
appears to have been a typical disk galaxy prior to the encounter. 

The ring shows electron densities consistent with typical values for star-forming 
HII regions. The eastern half of the ring exhibits a metallicity a factor of ~ 2 higher 
than the western half. The ionisation parameter, measured across the ring, roughly 
follows the previously observed trend with metallicity. 

Key words: galaxies - galaxies: individual: Arp 147- instrumentatiomspectrographs. 



1 INTRODUCTION 



1.1 Ring Galaxies and Arp 147 



Arp 147 is a ring galaxy with an early-type compan- 
ion from the Arp catalogue of peculiar galaxies flArp[ 1966). 
Taking H = 73 ± 3kms _1 Mpc _1 (jSpergel et al.|2Q07|) and 



Ring galaxies are a rare type of galaxy commonly observed 
with a nearby companion. Examples of such objects include 
the C artwheel galaxy (|Fosbury &; Hawarden||1977 Toomre 



given the measured systemic velocity of 9656 =b 18kms -1 
( |de Vaucouleurs et al.|p~991| ) the system is at a distance of 
132 Mpc implying a scale of 0.64 kpc per arcsecond. Some 



1978| |Davies fe Morton||1982| |Marcum et al.||1992), II Hz basic information for Arp 147 is shown in Table [T] 



4 flLynds fe Toomre|1976| |Theys fe Spiegel) 1976| |Appleton 



fe Marston|1997 MRomano et al.||2008), VII Zw 466 (|Theys 
k, Spiegel||1976| [Appleton Marston||1997| |Romano et al. 
2008) and AM 064-741 ( |Few et al.|1982p toiiame a few that 



Most ring galaxies are thought to be produced by a col- 
lision between an intruder galaxy and a normal disk galaxy, 



have been studied extensively. 
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as first put forward by |Lynds &; Toomre ( 1976| ) and |Theys k\ 
Spiegel (1977). In this scenario the intruder galaxy collides 



with a normal disk galaxy causing a gravitational distur- 
bance in the disk. As the intruder passes through the disk 
material is pulled towards the centre of the disk (i.e. the po- 
sition of the intruder) . The material then rebounds outwards 
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once the intruder has passed through the disk, causing an 
over density of material in the form of an outwards propa- 
gating ring. 

The basic parameters governing the geometry of the col- 
lision are the impact parameter and the angle of impact. The 
impact parameter describes how close the centres of the col- 
liding galaxies pass each other, measured as the percentage 
of the disk radius between the point of impact and the centre 
of the disk. The angle of impact measures how far from per- 
pendicular the collision is. Therefore in a purely symmetric 
collision, where the galaxies meet each other head on and 
perpendicularly, the impact parameter is and the impact 
angle is 0° 

The symmetric scenario was modelled very successfully 



Appleton & Marston (1997) found that ring galaxies 



by Mihos & Hernquist (1994) who examined an interaction 



between a disk of stars and gas in a dark matter halo and 
an intruding gasless system. Examining a head-on, perpen- 
dicular collision between these two objects they were able to 
recreate a system very much like the Cartwheel galaxy, with 
the nucleus of the original disk remaining in place. These au- 
thors also modelled the star formation triggered in the ring 
finding that the maximum star formation enhancement oc- 
curs late in the interaction, some time after the companion 
has passed though the ring. 

It was noted, however, by |Lynds &; Toomrej ( [1976) and 
Toomre ( 1978) that this highly symmetric collision can seem 



contrived. In fact these authors suggest that varying the an- 
gle of impact by up to 30° to normal can still yield ring-like 
shapes, also suggesting a modest impact parameter (up to 
15% of the disk radius) will likewise still produce ring-like 
structures. The variation of these parameters was studied 



in more detail by |Huang fc Stewart (1988) who modelled 



collisions with a combination of various impact angles and 
impact parameters. They found that variations of the impact 
angle of up to 30° show little deviation from a perpendic- 
ular impact with the nucleus of the disk not significantly 
dislodged. Larger impact angles and the introduction of a 
non-zero impact parameter, on the other hand, can cause 
the nucleus to become merged with or even offset from the 
resulting ring. These authors also conclude that an impact 
angle of 45° is most disruptive to the disk and that arc- like 
structures can still be formed for impact parameters up to 
about 50% of the disk length. 



In the context of Arp 147 Gerber et al. (1992) investi- 



gated an off-centre, perpendicular collision between systems 



similar to those used by Mihos & Hernquist (1994). They 
recreated a system like that of Arp 147, identifying the red- 
der region in the south-east of the ring as the nucleus and 
predicting it to be displaced relative to the plane of the 
ring. The IFS observations presented in this paper allow us 
to study the kinematics of the system in a spatially resolved 
manner and to examine the parameters of the collision that 
created the ring. 



It was predicted by Theys & Spiegel ( 1976) that the ring 
is a transient phenomenon on a time-scale of ~ 10 8 yrs. Since 
star formation is expected to be triggered by the expanding 
ring, there should be observational evidence of enhanced star 
formation in the ring and the star forming areas should have 
similar colours and ages. This makes ring galaxies very at- 
tractive to study since their star formation history should 
demonstrate a single spatially contiguous starburst around 
the ring. 



are generally blue and show a radial colour gradient in the 
sense that the colours become bluer in the outer part of the 
ring. This reinforces the idea that the star formation oc- 
curs in the ring and as the ring propagates outwards the 



stars age in its wake. Bransford et al. (1998) also showed 



that star-forming knots seen in most ring galaxies tend to 
have very similar colours, further suggesting that they have 
similar ages and were formed in the ring at the same time. 
Marston & Appleton ( 1995| ) found that ring galaxies have 
similar Ha luminosities to starburst galaxies, confirming the 
results of Appleton &; Struck-Marcell ( 1987) who found that 
ring galaxies have far infrared luminosities brighter than 
normal galaxies and comparable starburst galaxies. These 
authors also found that the star formation is almost exclu- 
sively found in the ring itself and not inside it. 

Some kinematic studies of ring galaxies have also been 



carried out. The pioneering study by Theys & Spiegel ( 1976) 



presented kinematics for two galaxies, VII Zw 466 and Arp 
147. For the latter system they were able to derive only 
tentative velocities which none the less were shown to be 
suggestive of expansion across the ring. For VII Zw 466 the 
authors were able to derive expansion and rotation veloci- 
ties for the ring by using data from many slit spectra taken 
at different positions across the ring. Fosbury & Hawarden 



(1977) carried out similar work for the Cartwheel galaxy, de- 
riving rotation and expansion velocities for the ring. Other 



studies include [Few et a~ (1982) who looked at slit spectra 
of AM 064-741 and |Taylor fe Atherton| fll984| who stud- 
ied the Vela ring galaxy using a Fabry-Perot interferometer. 



Appleton & Struck-Marcell ( 1996 ) presented some sample 



velocity maps from a ring galaxy simulation, which could be 
compared to observed two-dimensional kinematic maps. 

It is very clear that three-dimensional IFS data sets 
yielding two-dimensional kinematic maps of ring galaxies 
would be a rich addition to the field. The work presented in 
this paper provides an example of this. 

Some recent studies in various wavebands have fo cussed 
on Arp 147. Rappaport et ah] ( |2010| looked at muliti- 
wavelength images of the system, including infra-red and 
X-ray data. They found a star formation rate in the ring con- 
sistent with 7M0yr -1 , and derived an interaction timescale 
of ~ 40Myrs. They also located several ultra-luminous X- 
ray sources associated with the star-forming ring. 



Romano et al. ( 2008 ) also looked at Arp 147, deriving a 
star-formation rate of 8.6M©yr _1 from their measured far 
infra-red flux. From their measurements of the combined Ha 
and [Nil] flux they derive a star formation rate and equiv- 
alent width slightly lower than ours. Apple ton &; Struck- 1 
|Marcell ( 1987) derived a far infrared luminosity for Arp 147, 
corresponding to a SFR of 3.7M0yr _1 according to the re- 
lation oflKennicutt ( 1998 ). The star formation rates derived 



in these studies are broadly consistent with ours, given the 
scatter in the various relations used to derive the rates. 

The Arp 147 system, as well as being a fascinating ob- 
ject in its own right, is also an excellent example of an un- 
filled ring with a displaced nucleus. Several systems like this 
are known and an insight into their nature may result in 
a better understanding of the formation of these intriguing 
systems. 
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Arp 147 


RA (h:m:s J2000) 


03:11:18.9 


Dec (°:':" J2000) 


01:18:53 


Redshift 


0.032 


Systemic Velocity (kms -1 ) 


9656 ± 18 



Table 1. Arp 147 Basic Information 



1.2 The SWIFT Integral Field Spectrograph 
(IFS) 

The SWIFT Integral Field Spectrograph (IFS) is an instru- 
ment built by a dedicated team at the University of Ox- 
ford flThatte et aT1|2006| ). SWIFT is an AO-assisted, slicer- 
based IFS operating in the I & z bands. The instrument is 
mounted at the Cassegrain focus of the 200-inch Hale tele- 
scope at Palomar Observatory in California, where it was 
installed and commissioned in October 2008. SWIFT oper- 
ates with the Palomar Adaptive Optics system (PAL AO) 
in three modes - using the laser guide star (LGS-AO), us- 
ing a natural guide star (NGS-AO), or in static correction 
mode. The latter puts a fixed pattern on the deformable 
mirror to remove static aberrations in the telescope and AO 
system and provides no further correction. SWIFT has an 
on-the-fly selectable spatial pixel (spaxel) scale consisting of 
coarse, intermediate and fine spaxels (C//235, (//160, O'/OSO 
respectively). The instrument field-of-view is 44x89 spax- 
els in size yielding an on-sky field-of-view of 10 / /3x20 / /9, 
7 / /0xl4 / /2, dUbxf/l in the three spaxel scales. The instru- 
ment records ~ 4000 spectra in each exposure. There are two 
selectable wavelength ranges of 650nm-1050nm or 750nm- 
1050nm which are selected by changing a dichroic element in 
the PAL AO system. The instrument has a twin spectrograph 
design such that the field of view is split into two pseudo-slits 
at the image sheer. Each pseudo-slit then passes through 
a separate but identical spectrograph and is imaged onto 
one of two CCD detectors. SWIFT uses two 4kx2k, 250/xm 
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thick, fully depleted, CCD arrays developed by LBNL (Hol- 
2006). These devices have very high quantum 



efficiency up to ~ 1/im, even when compared to deep deple- 
tion devices ( [Burke et al~1|2004| ). 



2 THE DATA 

2.1 SWIFT Observations 

Arp 147 was observed with the SWIFT IFS on 12th January 
2009. The instrument was operated in coarse spaxel mode, 
with 0"235 spaxels, giving a field of view of 10 / /3x20 / /9. The 
wavelength range was 650nm-1050nm covering the region of 
the Ha; [Nil] AA 6548, 6583 A; [SII] AA 6716, 6731 A; and 
[SIII] AA 9069, 9532 A emission lines at the redshift of this 
object. 

Three pointings, labeled A, B, and C in Figure [I] were 
used to form a mosaic covering the entire system. The cen- 
tral coordinates of the exposures are given in Table [2] and 
the field of view is aligned with a PA of 90° (i.e. with the 
long axis of the field of view aligned N-S). To point to the 
object we first centered the SWIFT field of view on a bright 
star and offset to the desired position. The coordinates of 




Figure 1. A HST/WFPC2 image of the Arpl47 system with the 
observed SWIFT pointings labeled A, B, and C. The image is a 
false-colour composite of images taken in F450W, F606W, and 
F814W filters. Credit: NASA, ESA, M.Livio (STScI). 



Position 


RA 


Dec 


Tot. Exp. Time (s) 




(h:m:s J2000) 


J2000) 




A 


03:11:18.62 


01:18:56.8 


600 


B 


03:11:18.08 


01:18:57.9 


600 


C 


03:11:19.54 


01:18:48.1 


600 


Star 


03:11:24.81 


01:16:13.7 


N/A 



Table 2. Coordinates and Exposure Times 

the offset star are also given in Table [2] Two 300s- long ex- 
posures were taken at each pointing with an offset sky expo- 
sure taken between each object frame, giving an on-of f-on 
trio of frames for each pointing. A photometric standard 
star (%2 Cet) was observed after the science data was taken. 
Since this target was part of a bad weather back-up program 
these data were taken when conditions were less than opti- 
mal. During the observations the average seeing was mea- 
sured from SWIFT exposures of the bright offset star to be 
~ 2 / . / 5. There was also some cirrus cloud interfering with the 
observations, leading to variations in the sky background 
between observations. This meant that extra care had to be 
taken when reducing the data. The reduction procedure is 
described in Section 03 

2.2 HST/WFPC2 Observations 

Hubble Space Telescope (HST) observations of Arp 147 are 
available in the Hubble Archives. The system was observed 
using the Wide Field Planetary Camera 2 (WFPC2) in Oc- 
tober 2008, programme number 11902 (P.I. Mario Livio, 
STScI). Images were taken with the F450W, F606W and 
F814W filters, which are similar to the B, V and I bands 
respectively. The data were retrieved from the archive hav- 
ing already been reduced using the STScI pipeline. We then 
flux- calibrate the images using the PHOTFLAM keyword 
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determined by the pipeline and stitch the exposures together 
using the IRAF task wmosaic from the stsdas package. 



3 DATA REDUCTION 

3.1 SWIFT integral field data 

We reduce the SWIFT data in several steps using different 
packages, including the dedicated SWIFT data reduction 
pipeline (R. Houghton, private communication). Since the 
instrument has two twin spectrographs (as mentioned in Sec- 
tion 



1.2 ), each with a separate detector, there are two frames 
to be processed for each SWIFT exposure taken. These cor- 
respond to one half of the field of view, split along the short 
axis, and are labeled master and slave. As far as possible we 
process the master and slave frames separately, though we 
stitch them together to recover the entire SWIFT field of 
view at the end of the reduction procedure. 

The basic reduction procedure is the same for all sci- 
ence frames (consisting of the object, sky and standard star 
frames). We prepare the files using a pipeline routine which 
measures and removes the bias level, trims the overscan from 
the frame, and flat-fields the data. At this stage we apply 
the L. A. COSMIC routine ( |van Dokkum|2001| ) to identify and 
flag cosmic rays in the data. The identified cosmic rays are 
replaced with NaN (i.e. not a number) values. 

We find two wavelength solutions (one each for mas- 
ter and slave frames) using Ar and Ne arc lamp expo- 
sures. We apply the resulting wavemaps to the prepared 
two-dimensional science frames, yielding a set of master and 
slave data cubes. 

We then re-normalise the data cubes along the spec- 
tral axis by dividing through with a low-order polynomial 
representing the flat-field lamp spectrum. We flux calibrate 
the cubes using the standard star (%2 Cet) observed on the 
same night. At this stage we also apply corrections for the 
differing throughputs of the two spectrograph channels. 

Sky subtraction is complicated by the presence of vari- 
able cirrus clouds during the observations, causing large vari- 
ations in sky brightness on short time scales. We can anal- 
yse emission lines that do not require accurate continuum 
strengths by fitting a low-order polynomial to the galaxy 
continuum (including any residual sky background). This 
works well for lines uncontaminated by night sky emission, 
such as the Ha and [Nil] lines. We extract spatially inte- 
grated spectra in this way and two are shown in Figure [2] 

The [SII] emission line at 6731 A is contaminated by 
a sky emission line over a substantial fraction of the ring 
galaxy. To correctly subtract the sky emission, we fit six sky 
lines (including the contaminating one) in the near vicinity 
of the [SII] doublet with Gaussian line profiles. This is done 
for the sky spectrum, and the relative strengths, wavelength 
ratios and line widths derived from the fit are determined. 
We then hold these relative parameters fixed and fit to the 
same lines (allowing the overall strength to vary), exclud- 
ing the line which contaminates the [SII] emission, in each 
object spectrum. We use the five other sky lines and the rel- 
ative parameters already determined from the fit to the sky 
spectrum to determine the fit parameters for the contam- 
inating sky line. We then subtract this final fit to the sky 
lines, including the reconstructed contaminating line, from 
the object spectrum. The results are shown in Figure [3] 
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Figure 2. Spatially integrated spectra of the ring galaxy (a) and 
the early-type companion (b) showing the Ha and [Nil] lines. 
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a) SII lines with sky 
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Figure 3. The [SII] sky subtraction procedure. Panel a shows the 
raw spectrum covering the region containing the [SII] emission 
lines and sky lines, with the Gaussian fits to the sky lines super- 
imposed. Panel b shows the result of subtracting this fit from the 
raw spectrum. 



3.2 Colours and equivalent widths 

We use the HST data to measure spatially resolved colours 
for the ring galaxy and the early type galaxy. Multiple expo- 
sures are available in each of F450W, F606W, and F814W 
filters. We combine all exposures in a single filter after re- 
alignment and rejection of cosmic rays. The resulting data 
sets are flux calibrated using the latest PH0TFLAM values, 
and registered using world co-ordinate (WCS) information. 
We then made a mask frame based on a minimum surface 
brightness cut-off in the F606W image. An F450W-F814W 
(approximately B— I) colour map created using the masked, 
re-aligned, flux- calibrated data sets is shown in Figure[4](the 
error between B— I colour and observed F450W— F814W 
colour depends on spectral type, but is typically about 0.02 
mag). We also determined the mean colour in each of four 
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Arpl47 Velocity Map 



0.0 0.5 1.0 1.5 2.0 2.5 

Figure 4. F450W-F814W colour for the Arp 147 system. The 
colour values have been corrected to correspond directly with B— I 
colours in the Vega system. The western part of the ring galaxy 
has significantly bluer colours, while the reddest colours are in 
the SE region, possibly corresponding to the nuclear bulge of the 
original disk galaxy from which the ring was formed. 



quadrants of the ring galaxy (defined in Section 4.2), the 
measured values are tabulated in Table [5] Typical errors on 
the colours are less than 0.05 mag. 

To measure the Ha equivalent width it is necessary to 
have some estimate of the galaxy continuum at wavelengths 
close to the emission lines. We measure the continuum level 
from integrated HST photometry in the F606W filter, which 
has a pivot wavelength of 5919 A and a bandwidth of 
1579 A. Given the worse-than- average seeing conditions dur- 
ing our observations, we restricted the spatial resolution of 
our equivalent width measurements to a single value for each 



quadrant of the ring galaxy (see Section 4.2 and Figure [9j). 
We mask the HST image to exclude any sky, and then sum 
the flux in each quadrant. We then subtract from each quad- 
rant the measured line fluxes in the Ha and both [Nil] lines 
in that quadrant, to correctly account for the line contami- 
nation. This allows us to compute an estimated continuum 
level in each quadrant, in units of ergs -1 A -1 cm -2 , and the 
equivalent width for each quadrant. The same procedure is 
followed to derive a single integrated value for the early type 
galaxy. 



4 RESULTS 

4.1 Collision Geometry 

We fit an ellipse to the image of the ring galaxy in the F814W 
band and use the fitted major and minor axes of the ellipse 
to find the inclination of the ring using Equation [T| ( |Fosbury| 
&; Hawarden||1977 ), which assumes an intrinsic circular ge- 
ometry. 
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Figure 5. The velocity map of the Arpl47 system. This is derived 
by fitting the Ha emission line in the ring galaxy with a Gaus- 
sian and thus calculating the corresponding line-of-sight velocity 
shift. The same procedure is used for the early- type galaxy but 
instead fitting the comparatively stronger [NII]A6583 line. The 
kinematic PA of the ring, measured as described in Section |4~T] 
is superimposed. 
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We find that the inclination of the ring is 25.3 =b 1.5° 
and the radius of the ring is 5.8 ±0.1 kpc. The photometric 
position angle (PA) of the ring is 33 =b 4° where the errors 
on these values are derived from the fitted ellipse. 

We obtain a velocity map of the Arp 147 system using 
the observed emission lines. For the ring galaxy we isolate 
and fit the Ha line with a Gaussian profile in each spaxel 
where the line has a S/N^3.5. We then use the centroid 
of the fit to find the line-of-sight velocity. The procedure is 
the same for the early-type galaxy, but we use the stronger 
[Nil] emission line. Figure [5]shows the measured line-of-sight 
velocity map for the entire system. 

For the early-type galaxy we do not expect a signifi- 
cant line-of-sight velocity gradient, since the emission used 
to calculate the velocity is not stellar in origin (see Section 
4.2). Since we lacked stellar absorption lines to determine 
the velocity this emission line velocity is included in the 
map merely for completeness. 

The observed line-of-sight velocity gradient across the 
ring galaxy is due to rotation and expansion of the ring. 
These two components must be disentangled from each 
other. As shown in Figure [6] we do this by measuring the 
photometric and kinematic position angles (PAs) of the ring. 
Assuming an intrinsic circular geometry, if the ring is purely 
rotating we expect to measure the kinematic and photomet- 
ric PAs to be the same. This is illustrated in Figures |5Jl and 
^Jd. If the ring is purely expanding we expect to measure the 
kinematic and photometric PAs to be perpendicular to each 
other. This is illustrated in Figures [5J: and[6]l. Therefore if 
we calculate the angle between the measured kinematic and 
photometric PAs of the ring galaxy it is possible to calcu- 
late how much of the measured line-of-sight velocity is in 
rotation, and how much is in expansion. We measure the 
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Figure 6. The ring galaxy geometry. Panels a and b show that a 
purely rotating ring will exhibit a kinematic major axis coincident 
with its photometric major axis. Panels c and d show that a purely 
expanding ring will exhibit a kinematic major axis perpendicular 
to its photometric major axis. 



kinematic PA of the ring galaxy as 70.5 =b 0.5° using an IDL 
routine (by M. Cappellari, private communication) which 
uses the method described in Appendix C of Krajn ovic et al.| 
( 2006 ) . The kinematic PA is shown superimposed on the ve- 



locity map in Figure [5] 

To measure the component of the observed velocity due 
to expansion we plot each measured velocity point as a func- 
tion of the azimuthal angle of the ring, 9. We then fit the re- 
sulting graph with Equation |2|( |Fosbury &; Hawarden|1977 ), 
where V sys t is the systemic velocity of the ring, V ro t the com- 
ponent of the observed velocity due to rotation and V exp 
component due to expansion of the ring. 



the 



V = Vsyst + sin(i)(Vex P sin(6>) + V ro tcos(6>)) 



(2) 

The deprojected edge-to-edge expansion of the ring, 
Vexp, is 225 ± 8 km s -1 , giving a radial expansion velocity of 
113 ± 4kms -1 . Given the radius of the ring (5.8 =b 0.1 kpc) 
and assuming a uniform expansion velocity from a point at 
the centre of the ring, this implies that the collision occurred 
a maximum of 50 Myrs ago. Comparatively, there is little ro- 
tational velocity, with a value for V ro t of 47 =b 8kms _1 . 

We also measure the systemic velocities of the ring and 
the early-type companion by fitting emission lines in inte- 
grated spectra of each object. The difference between the 
line-of-sight velocity of the two objects is 151 ± 7kms _1 . 
Given the upper limit on the collision timescale found above 
we use the age of the minimum age of the interaction- 
triggered young population, taken to be 1.1 Myrs (see Sec- 
and Table [5]) to calculate 



tion 



4.2.1 



lower limit for the 

interaction timescale. 

Since the star formation was more than likely delayed 
from the time of impact ( Mihos &; Hernquist| 1994) we arbi- 
trarily assume it took a minimum of 10 Myrs for star for- 
mation to be triggered (Rappaport et al. 2010) and use 



Property 


Value 


Ring Inclination 


25.3 ± 1.5° 


Ring Photometric PA 


33 ±4° 


Ring Kinematic PA 


70.5 ±0.5° 


Ring Radius 


5.8±0.1kpc 


Ring Expansion Velocity 


225 ± 8 km s- 1 


On Sky Separation a 


12.8 kpc 


Maximum Interaction Timescale 6 


50 Myrs 


Impact Angle Range 


33° - 54° 



Table 3. Kinematic Properties of the Arp 147 System. a On sky 
separation between the ring and companion. b Upper limit on 
interaction timescale as described in the text. 



11.1 Myrs as a lower limit on the interaction timescale. This 
implies that the early-type companion has travelled between 
1.7 =b 0.1 kpc and 7.8 =b 0.5 kpc along the line of sight since 
the collision occurred. The distance between the centres of 
the galaxies in the plane of the sky is 12.8 =b 0.1 kpc which 
implies that the collision occured at an impact angle of be- 
tween 33° and 54°. A 0° impact angle implies a head-on 
collision. 



4.2 Star Formation 

The SWIFT Ha and [Nil] A6583 linemaps for the Arpl47 
system are shown in Figures [T^l and [7J) respectively. The 
maps are continuum subtracted and the wavelength ranges 
used are centred on the lines in question (with different 
ranges for each object accounting for differences in line-of- 
sight velocity). The ring galaxy is obviously strongly star- 
forming. The source of ionisation in the early- type compan- 
ion is less obvious, however. Typically emission line objects 
are classified by comparison of their line ratios through the 
use of so-called BPT diagrams Baldwin et al. (1981). Fig- 
ure [8] shows the parameter space of a typical BPT diagram 
which compares the log([NII]/Ha) and log([OIII]/H/3) ra- 
tios. Our spectra do not enough lines to place the objects in 
the Arpl47 system onto a BPT diagram, but it is possible 
to distinguish the mechanisms behind the line emission in 
these objects none the less. 

|Kewley et al.| ( |2001| used stellar evolution models to 
derive a theoretical limit for starburst galaxies on the BPT 
diagram. This line is shown as the green curve in Figure [8] 
and represents an extreme upper limit for starburst galax- 
Kauffmann et al. ( 2003[ ) revised this downwards to sep- 



arate AGN and star-forming galaxies in an empirical way. 
Their line is shown in Figure [8] as a dashed cyan line. Ob- 
jects between the two lines are considered to be "composite" 
objects whose ionising spectrum arises due to a combina- 
tion of star formation and AGN. For the Arp 147 galaxies 
we can measure only one parameter of the BPT diagram - 
log([NII]/Ha). The ring galaxy is shown as a blue line along 
its possible log([OIII]/H/3) values. It is clear that this ob- 
ject is strongly star- forming (not that there was any doubt). 
The early-type companion is shown as a red line covering 
its possible log([OIII]/H/3) values. This galaxy clearly falls 
in the LINER/Seyfert part of the BPT diagram. 

We extract spatially integrated spectra for each of the 
two galaxies in the system. We also study the ring in more 
detail to see if any of the physical properties vary over its 
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Figure 7. a) The Ha linemap of the Arpl47 system, b) The 
[Nil] A6583 linemap of the Arpl47 system. Both maps have been 
continuum subtracted. 



extent - we extract four spectra from quadrants across the 
ring (see Figure [9| . We chose this approach to ensure the 
spatially integrated spectra had enough S/N to measure the 
weaker emission lines. The quadrants are not equally sized 
and were chosen such that quadrant 1 covered the redder 
region to the south-east of the ring. It is possible that this 
region corresponds to the nuclear bulge of the original disk 
galaxy from which the ring was formed. 

We measure the integrated Ha flux for the ring galaxy 
and use this to derive the star formation rate according to 
Equation [3] ( Kennicutt||1998t . We measure a star formation 
rate of 4.68±0.67 M /yr, as shown in Table 



SFR(M©yr~ 



7.9 x 10" 42 L(Ha)(ergs" 



(3) 



4-2.1 Star Formaton history 

We measure the Ha equivalent width for each quadrant of 



the ring galaxy (see Section 3.2), using spatially integrated 
spectra from each region. The results are presented in Ta- 
ble [5] Using the measured equivalent widths and colours, we 
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Figure 8. A BPT diagnostic diagram showing the 'extreme star- 
burst line' (green ref) and the empirical line below which all galax- 
ies are starbusts (dashed cyan ref). The red line vertical line rep- 
resents the early-type galaxy, consistent with a non-star-forming 
ionisation mechanism (either a LINER or a Seyfert galaxy) and 
the blue line represents the star-forming ring galaxy 




Figure 9. The ring galaxy was split into the four quadrants 
shown, labelled 1, 2, 3 and 4. Quadrant 1 was chosen to include 
the redder area in the south-east corner, thought to correspond 
to the nuclear bulge of the original disk galaxy. A spatially inte- 
grated spectrum was extracted for each quadrant. 
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Figure 10. Plot showing the range of allowable ages for the recent 
star formation in three of the four quadrants of the ring galaxy. 
There are two curves for each quadrant (black for Q2, red for Q3 
and blue for Q4), representing the extremes of the allowed val- 
ues for the Ha EW as a function of age of the young (encounter 
triggered) stellar population. These extremes take into account 
instantaneous and continuous star formation histories, and all al- 
lowable ages for the older (pre-encounter) stellar population in 
the disk galaxy. The observed value of the Ha EW in each quad- 
rant defines an age range for the young stars, these are shown by 
horizontal bars for each quadrant. 



can constrain the age range of the recent star formation ac- 
tivity occurring in the ring galaxy, by comparing with stel- 
lar population synthesis models. In particular, we use the 



Starburst99 (Leitherer et al. 1999) code, with a standard 



Salpeter ( 1955 ) initial mass function (IMF) with upper and 
lower mass cutoffs of 0.1 M© and 100 respectively, and 
solar metallicity. We consider both instantaneous starbursts 
and continuous star formation as the two limiting cases for 
the young (encounter triggered) stellar population. This is 
consistent with |Mihos &; Hernquist (1994), who predict an 
enhancement in star formation in the ring of a factor of 2 — 3. 

Our analysis to derive the age range of the recent star 
formation is best understood as a step-by-step process. 

Step 0: Assumptions It is reasonable to expect that the 
observed light from the ring galaxy consists of two compo- 
nents, one stemming from an old population (native stellar 
population of the disk galaxy prior to the encounter) , and a 
second younger population, triggered by the encounter. We 
assume that the old population native to the disk galaxy has 
the same age in all quadrants. Also, we assume negligible ex- 
tinction in all quadrants of the ring. The latter assumption 
is based on the visual appearance of the ring (see Figure [I]) 
and the colour map (Figure [4]), the low value of the inclina- 
tion, and is partly motivated by the lack of sufficient data 
to derive a reliable extinction estimate. The reddest quad- 
rant (Ql of Figure [5]) which appears to contain the nuclear 
bulge of the original disk galaxy may contain a stellar pop- 
ulation of a different age (bulge component), and is treated 
separately. 

Step 1: Constraining the age range of the old stellar pop- 
ulation The reddest B— I colour we observe in a non-nuclear 
quadrant (Q2) is 1.6, so the age of the old component of 
the disk galaxy is restricted to correspond to B— I colour 



redder than 1.6 (the colour gets bluer when we add a young 
component). That gives a lower bound on the age of the 
old component of 1.3 Gyr, the upper bound is the age of 
the Universe, or 13 Gyr. (We note that the B— I colour of 
the neighbouring early type galaxy is consistent with a stel- 
lar population of 13 Gyr age). Knowing the age (and hence 
colour) of the old component, we can determine the fraction 
of light from young stars. The observed Ha equivalent width 
(EW) in each quadrant places a further constraint on the 
old component colour, as too small a fraction of young stars 
cannot simultaneously account for the observed Ha EW and 
B— I colour. This increases the lower limit on the old compo- 
nent to 1.8 Gyr, if the young stars form in an instantaneous 
starburst, or 1.5 Gyr if they form continuously. 

Step 2: Fraction of light from young stars For a fixed old 
component age, we can compute the young population frac- 
tion required to produce the observed colour in each quad- 
rant, as a function of age of the young component. Two 
curves are computed, one for the young stars forming in a 
single burst (instantaneous) or forming continuously. Vary- 
ing the old component age over the allowed range gives us 
two families of curves showing the fraction of light from 
young stars, one for instantaneous and the other for con- 
tinuous star formation. 

Step 3: Computing the Ha EW v/s young component 
age Each curve of light fraction v/s age of the young com- 
ponent (for a fixed old population age) corresponds to a 
curve of observed Ha EW v/s age of young component. Fig- 
ure ^] show the most extreme curves for three quadrants 
(excluding the nuclear quadrant) - black for Q2, red for Q3, 
and blue for Q4. 

Step 4 : Constrain age and light fraction using Ha EW 
For an observed Ha EW, the region between the most ex- 
treme curves corresponds to the allowed age range for the 
young component in each quadrant, this is shown by hori- 
zontal bars in Figure [l0| We get an age range, rather than 
a fixed value, due to the uncertainty in knowing the age of 
the old component, and the recent star formation history 
(i.e. instantaneous or continuous). For each quadrant, we 
can derive an allowed age range and the corresponding frac- 
tion of light from the young star formation triggered by the 
encounter. These are presented in Table [5] The age range for 
the young stars is modest, ranging between 1.1 and 70 Myr, 
independent of star formation history. We use the mass-to- 
light ratio for each component provided by the Starburst99 
models to also derive a mass fraction in young stars. 

It is interesting to note that although the young stars 
produce as much as a third of the V band light in the ring 
galaxy Q4, the mass fraction is negligible at ^1% in all cases, 
irrespective of which recent star formation mode we adopt 
- a robust result. We also note that all quadrants, including 
those dominated by blue light from young stars require a 
two component stellar population to fit the measured values 
for the colour and equivalent width simultaneously. 

In quadrant 1 of the ring galaxy, the observed colour is 
too red (B— I = 2.09) to allow a significant age range for the 
old population. Indeed, the visual appearance (see Figure 
|9| suggests that some reddening may be due to dust. For 
consistency with the other measurements, we compute an 
age range for Ql, assuming that the old population is 13 Gyr 
old. The results are shown in Table [5] The corresponding 
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Object 


L(Ha) (erg/s) 


SFR (M /yr) 


Ring 


5.92(±0.85) x 10 41 


4.68 ±0.67 



Table 4. The Ha flux and star formation rate for the ring galaxy. 

light fraction for the young component is very small, and 
the mass fraction is entirely negligible, as expected. 

We do not do this computation for the early type galaxy, 
as the observed Ha emission has a different origin, as dis- 
cussed in Section l4~2l above. 



4-2.2 Total stellar mass 

Using the results of the previous section for the light frac- 
tion from the old component in each quadrant of the ring 
galaxy, and using the total observed F606W continuum flux 
from the ring galaxy, we compute the total stellar mass of 
the ring galaxy to be in the range 6.1 — 10.1 x 10 10 M©, the 
uncertainty stemming from the age of the old stellar compo- 
nent that dominates the mass. Thus, it appears as though 
the ring galaxy is a typical disk galaxy, similar to the Milky 
Way. The "nuclear bulge" (Ql) comprises 30-50% of the to- 
tal mass. 



4.3 Line Ratios and Derived Properties. 

The N2 ratio (N2 = log([NII]/Ha)) is a relativ ely reliable 
empirical indicator of oxygen over- abundance ( |Pettini "fc 



|Pagel|2004| . This relation is justifiable because at the metal- 
licity range studied here oxygen can be assumed to be a 
product of primary nucleosynthesis and nitrogen a product 
of secondary nucleosynthesis. Therefore the strength of ni- 
trogen emission gives an indication of the over- abundance 



of oxygen in the star- forming region ( Storchi-Berg mann| 
et al.|1994 ). We calculate oxygen overabundance for the ring 
galaxy and for each quadrant of the ring using Equation [4] 
where the solar value is 12+log(O/H)=8.76±0.07 flCaffau| 
et al.|2008 ). The results are presented in Table [6] Of note is 
the fact that the eastern half of the ring shows a roughly so- 
lar oxygen abundance, whereas the western half of the ring 
shows an oxygen abundance a factor of two lower. 



12 + log(0/H) = 8.9 + 0.57 • log([NII]/Ha) 



(4) 



Since the [SII] AA 6716, 6731 emission lines have simi- 
lar excitation energy but differing collisional de-excitation 
rates, the observed emission strengths of these two lines 
are dependent on the electron density in the emit- 
ting region ( |Osterbrock| [1989). We use the ratio S2 = 
[SII](A6716)/[SII](A6731) to calculate the electron density, 
n e , in each quadrant of the ring galaxy and in the early- 
type companion. We use the standard iraf nebular package 
to find the values which are shown in Table [6] All quadrants 
of the ring galaxy have typical electron densities expected 
for star-forming HII regions ( Kennicutt et al.|1989 ). 

The ionisation parameter, u, is the ratio of the ionising 
photon density to the particle density in the HII region and 
is one of the parameters defining the emission line spectra 
of HII regions. It can be found using the measured ratio 
[SII] (6717 + 6731 )/[SIII](9069 + 9532) and Equation [s| (see 
piazet al.|p99H ). The value of log(u) was calculated for the 



ring galaxy and each quadrant of the ring. The calculated 
values of log(u) are shown in Table [6] 

log(u) = -1.684 • log([SII]/[SIII]) - 2.986 (5) 
The ionisation parameter is expected to fall with in- 



creasing metallicity (Evans & Dopita 1985), and our mea- 



surements roughly follow this expected trend with oxygen 
abundance. 



5 CONCLUSIONS AND DISCUSSIONS 

The ring galaxy in the Arp 147 system is expanding with 
a radial velocity of 113 =b 4kms -1 . Assuming the expansion 
is a result of the impact which originally formed the ring, 
and assuming the expansion velocity remained constant and 
all material in the ring originated at a point in the centre 
of the ring, this implies an upper limit on the timescale for 
the collision of 50Myrs. The relative velocity between the 
early- type galaxy and the ring is 151 ± 7kms _1 . 



Mihos fc Hernquist] (|1994|) have created models of col- 



lisionally created ring galaxies by studying centred head-on 
collisions. Using the measured kinematics of Arp 147 it is 
shown that the impact which created Arp 147 probably oc- 
curred at an angle in the range 33° —54°. Evidence presented 
by |Gerber et al.| ( |1992[ ) implies that an off-centre collision 
better produces an empty ring with a dislodged nucleus, like 
that seen in Arp 147. Thus we speculate that an off-centre, 
angled collision created the Arp 147 system. The parameters 
we derive for the impact are shown in Tableland could be 
used to constrain future models. In their simulation Gerber 



et al. ( 1992| ) found that the nuclear region of the original 
disk was not only dislodged during the collision but was off- 
set from the plane of the resulting ring. They identified the 
dislodged nucleus with quadrant 1 of the actual ring as we 
do, but we see no evidence in our kinematics to suggest that 
the nucleus is indeed offset from the plane of the ring. 

We measure the equivalent width of the Ha lines in 
the ring galaxy, and the B— I colours from multi-band HST 
photometry. Comparing this with the results from the Star- 



burst99 stellar population synthesis models (Leitherer et al 



1999), we conclude that we require a two component stellar 
population throughout the ring galaxy, even in the bluest 
quadrant whose visual appearance is dominated by light 
from young stars. We constrain the age range of the recent 
star formation activity in the ring, finding that the young 
stellar population of the ring galaxy is between 1.1-70 Myr 
old, with the most extreme assumptions for star formation 
history (instantaneous and continuous). We find the corre- 
sponding fraction of light from the young stars to span the 
range 0.04 to 0.34, while the stellar mass formed during the 
recent star formation activity is negligibly small, below 1% 
and likely -0.1%. 

The models of |Mihos &; Hernqu ist (1994) give a good 
description of the evolution of star formation in the cre- 
ated ring in the aftermath of the collision. They find that 
some time after the encounter star formation tends to be 
dominated by the created ring. Their simulations show that 
the created ring moves outwards, sweeping up gas and trig- 
gering star formation. The inner parts of the ring are left 
devoid of star formation. The star formation in the ring is 
therefore not triggered immediately upon collision with the 
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Object 


EW(Ha) 


B-I 


iogio 


Age 


Light 


Mass 




(A) 


colour 


(years) 


fraction 


fraction 


Ring Ql 


-83 ±2.5 


2.09 


6.06 - 


6.27 


-0.04 




Ring Q2 


-76 ± 1.6 


1.60 


6.26 - 


7.84 


0.05 - 0.2 


<1% 


Ring Q3 


-151 ±3.0 


1.55 


6.04 - 


6.79 


0.06 - 0.13 


<0.1% 


Ring Q4 


-203 ±5.8 


1.21 


6.35 - 


7.27 


0.14 - 0.34 


<0.15% 


ETG 


-4.1 ±0.3 


2.27 











Table 5. The measured Ha equivalent widths, young starburst age ranges, light and mass fractions for the quadrants of the ring galaxy. 



Object 


12+log(0/H) 


S2 


n e 




log(u) 


Ring 


8.44 ±0.04 


1.16 ±0.32 


201cm- 


-3 


-3.1 ±0.1 


Quadl 


8.50 ±0.07 


1.16 ±0.15 


297cm- 


-3 


-2.9 ±0.1 


Quad2 


8.28 ±0.24 


1.25 ±0.21 


201cm- 


-3 


-3.0 ±0.1 


Quad3 


8.52 ±0.07 


1.20±0.58 a 


267cm- 


-3 


-3.1 ±0.1 


Quad4 


8.19 ±0.13 


1.25 ±0.20 


153cm- 


-3 


-2.8 ±0.1 



Table 6. The S2 parameter, electron density and ionisation parameter. a The [SII] A6731 line in this quadrant was particularly badly 
contaminated by sky so it was necessary to fix the width of the line to match that of the [SII] A6716, resulting in a large error on the fit 
to the line. 



companion, but with a certain delay, which they estimate to 
be ~ 10 Myrs. In the case of Arp 147, where the encounter 
was likely off-centre and not head-on, the colour map sug- 
gests an age gradient in the recent star formation, with the 
bluest colours in the west- southwest region, and progres- 
sively getting redder counterclockwise. This observation is 
consistent with the scenario that the "ring" galaxy is not 
a planar configuration, but is actually a "corkscrew" shape 
that appears planar when projected along the line-of-sight 
(M. Livio, priv. comm). Further detailed modelling will be 
required to confirm or deny this hypothesis. 

The redder colours in quadrant 1 of the ring could be 
due to dust reddening or they may indicate that this region 
corresponds to the nuclear bulge of the disk galaxy from 
which the ring was created. We compute a total stellar mass 
of 6.1 — 10.1 x 10 10 Mq for the ring galaxy, making it a rather 
typical disk galaxy prior to the encounter. We observe that 
between 30-50% of the stellar mass is in Ql, lending weight 
to the interpretation that this quadrant contains the bulge 
of the disk galaxy. 

We measure the oxygen over- abundance for the ring 
galaxy and in the four quadrants of the ring. The eastern half 
of the ring has an oxygen abundance a factor of ~ 2 higher 
than that of the western half of the ring. The electron den- 
sity and ionisation parameter values for the quadrants of the 
ring are consistent with those expected in star-forming HII 
regions. 
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